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KINETICS  OF  THE  DISPLACEMENT  REACTION  BETWEEN  Fe  AND  Cu20 


AT  TEMPERATURES  BETWEEN  800  and  1050°C 

by 

Gregory  J.  Yurek*,  Robert  A.  Rapp**  and  John  P.  Hirth** 

ABSTRACT 

The  products  of  the  displacement  reaction  between  solid  iron 
and  CUjO  in  the  temperature  range  of  800  to  1050#C  are  solid  copper, 
wustite  and  magnetite.  The  product  oxide  in  contact  with  the  iron 
phase  is  in  the  form  of  rods  embedded  in  a matrix  of  copper.  A thin 
layer  of  copper  separates  the  oxide  rods  from  the  Cu^O  phase.  The 
rates  of  growth  of  the  product  phases  are  controlled  by  solid-state 
counterdiffusion  of  oxygen  and  iron  through  the  product  phases, 
except  for  the  initial  stage  of  the  reaction.  A previously  developed 
kinetic  model  for  the  diffusion-controlled  growth  of  the  product 
phases  is  employed  to  calculate  the  rates  of  growth  of  each  of  the 
product  phases.  Theoretical  values  of  growth  rates  are  in  excellent 
agreement  with  experimental  values  for  temperatures  above  950°C. 

Below  950*C,  the  mechanism  for  the  growth  of  the  magnetite  phase 
changes,  which  yields  marked  differences  between  the  experimental 
and  theoretical  rates  of  growth.  The  mechanism  for  the  displacement 
reaction  below  950"C  has  not  been  determined.  ‘ 
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1.  INTRODUCTION 


A displacment  reaction  between  solid  iron  and  solid  CujO 
results  in  the  formation  of  an  aggregate  of  solid  product  phases 
between  the  reactant  phases  [1,2].  The  product  phases,  which  have 
been  identified  as  wustite,  magnetite  and  copper  [3],  are  arranged 
in  two  zones  between  the  reactant  phases  as  shown  schematically  in 
Fig.  1.  One  zone  is  located  adjacent  to  the  iron  and  consists  of 
rods  of  the  product  oxides  in  a matrix  of  copper.  The  other  zone 
is  a layer  of  copper  that  is  located  between  the  copper  + oxide 
zone  and  the  CUjO.  This  type  of  product  arrangement  has  been 
classified  as  a lamellar-aggregate  product  morphology  [2] , although 
"rod-aggregate"  product  morphology  is  more  descriptive. 

Growth  of  the  product  phases  of  the  Fe/CUjO  displacement 
reaction  requires  solid-state  diffusion  of  iron  cations  through 
the  oxide  rods  and  of  oxygen  atoms  through  the  product  copper  layer. 
Electons  can  migrate  from  the  Fe/FeO  interface  to  the  tips  of  the 
oxide  rods  through  the  adjacent  copper  phase  [see  Fig.  1.].  The 
oxide  rods  grow  at  the  boundary  between  the  two  product  zones  (the 
reaction  front)  as  a result  of  the  simultaneous  arrival  of  iron 
cations,  electrons  and  oxygen  atoms.  The  copper  layer  grows  at  the 
Cu/CUjO  interface  where  Cu20  dissociates  to  yield  metallic  copper 
and  oxygen  dissolved  in  the  copper. 

Equations  have  been  derived  to  calculate  the  equivalent  fluxes 

of  oxygen  atoms  and  of  iron  cations  at  the  reaction  front  [3].  The 

equivalent  fluxes  depend  on  two  independent  variables:  (.i)  p!l  , the 
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oxygen  potential  at  the  reaction  front  (see  Fig.  1)  and  (ii)  f, 
the  volume  fraction  of  oxide  in  the  copper  + oxide  product  zone. 
Two  equations  are  necessary  to  solve  for  the  two  variables,  and 
hence,  for  the  equivalent  fluxes  at  the  front.  One  equation 
is  obtained  from  a mass  balance;  i.e.  the  steady-state  equivalent 
fluxes  of  iron  and  oxygen  at  the  reaction  front  are  equal. 


A second 


equation  is  obtained  from  the  postulate  [3]  that 


the  two  variables  adjust  such  that  the  maximum  equivalent  fluxes 
steady  state,  thereby  maximizing  the  reaction  kinetics.  An 
analytical  solution  for  the  two  variables  has  not  been 
found;  however,  values  for  the  equivalent  fluxes  can  be  computed 

p 

as  a function  of  f or  P_  after  one  of  the  variables  is  eliminated 

°2 

4 p 

from  the  flux  equations.  Theoretical  values  of  f and  Pn  can  then 

u2 

be  used  to  calculate  the  rates  of  growth  of  the  various  product 
phases.  Predicted  reaction  rates  for  the  Fe/Cu20  displacement 
reaction  at  1000 °C  were  found  previously  to  be  in  good  agreement 
with  experimental  reaction  rates  (3].  The  purpose  of  the  present 
paper  is  to  extend  the  experimental  and  calculated  results  for 
the  Fe/Cu^O  displacement  reaction  to  temperatures  between  800  and 
1050*C. 


2.  THEORETICAL  ANALYSIS 

The  equations  and  data  that  are  necessary  to  calculate  the 
equivalent  fluxes  of  iron  cations  and  oxygen  atoms  at  the  reaction 
front  (Eqs.  (34)  and  (35)  of  [3])  are  presented  in  this  section 
with  some  modifications. 


2.1  Phases  Present 


The  upper  llsilt  of  temperature  for  applicability  of  the 
kinetic  model  for  the  Fe/CUjO  reaction  is  1065*C,  which  is  the 
eutectic  tempera tu«. . between  copper  and  Cu^O.  The  formation  of 
oxides  in  addition  to  wustite  in  the  product  zone  depends  on  the 
kinetically  established  value  of  P^  at  the  temperature  of 
interest.  The  oxygen  potential  in  the  product  zone  varies  from 
that  for  equilibrium  between  iron  and  wustite  to  that  for 
equilibrium  between  copper  and  Cu20.  According  to  the  oxygen 
potentials  for  equilibria  among  the  compounds  in  the  Cu-Fe-0 
system  [4,5,23],  wustite,  magnetite  and  delafossite  (CuFe02> 
are  possible  product  oxide  phases  in  the  temperature  range  of 
800  to  1050°C.  The  oxygen  partial  pressure  for  equilibrium  between 
copper,  magnetite  and  delafossite  is  sufficiently  close  to  the 
partial  pressure  for  equilibrium  between  copper  and  cuprous  oxide 

p 

that  delafossite  is  not  a likely  product  phase;  i.e.  as  Pn 

u2 

approaches  P (Cu/Cu.O) , the  diffusion-controlled  growth  of  the 
°2 

product  copper  layer  must  approach  zero.  Because  delafossite 
was  not  observed  experimentally  it  is  not  considered  further  in 
the  analysis. 


2.2  Parabolic  Rate  Constants  for  Oxides 

Calculation  of  the  equivalent  fluxes  of  iron  cations  and 
oxygen  atoms  at  the  reaction  front  requires  an  evaluation  of  the 
parabolic  rate  constants  for  the  diffusion-controlled  growth  of 
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wustite,  kptFe^O),  and  of  magnetite,  kptFe^O).  The  Fe/O  ratios 
for  these  oxides,  v and  6,  are  taken  as  average  values  of  the  Fe/O 
ratios  over  their  ranges  of  stability.  The  Fe/O  ratio  for  magne- 
tite is  approximately  0.75  for  the  entire  temperature  range  of 
interest.  The  Fe/O  ratio  for  wustite  exhibits  a considerable  vari- 
ation with  oxygen  potential  at  each  temperature  in  the  range  of 
800  - 1050 °C.  The  average  values  of  v given  in  Table  I were  cal- 
culated from  the  data  for  wustite  compiled  by  Giddings  and  Gordon 
(51.  The  use  of  an  average  Fe/O  ratio  for  wustite  is  justified 
on  the  basis  of  a recent  paper  (6]. 

The  parabolic  rate  constant  for  the  diffusion-limited  growth 
of  an  oxide  Me&0  on  a metal  Me  in  an  oxygen-containing  gas  is 
given  by 
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k (Me  0) 
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,P°2  <0L%e)dl"P0. 


where  D*  is  the  tracer  diffusion  coefficient  for  Me  cations  in 
Me 

Me  0,  g is  the  correlation  factor  for  cation  diffusion  in  Me  0 

01  M6  0l 

and  Pi  and  P"  are  the  oxygen  partial  pressures  at  the  metal/oxide 
2 °2 

and  oxide  gas  interfaces,  respectively.  Equation  (1)  is  a 
slightly  modified  form  (7)  of  the  equation  derived  by  Wagner 
[8,9]  for  the  parabolic  rate  constant.  The  major  assumptions 
involved  in  the  derivation  of  Eq.  (1)  are  that  (i)  the  variation 
in  stoichiometry  for  Mea<5  is  small,  (ii)  MeQ0  is  an  electronic 
conductor  and  (iii)  the  mobility  of  the  cations  is  much  greater 
than  the  mobility  of  anions  in  Mea0. 
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These  assumptions  apply  for  magnetite  with  a - 0.75. 
Assumptions  (ii)  and  (iii)  apply  for  Fe^O.  Assumption  (i)  is 
not  strictly  applicable  for  FeyO,  but  the  use  of  the  average 
values  v at  the  temperatures  of  interest  is  justified  on  the 
basis  of  the  results  in  ref.  (6]. 

When  both  wustite  and  magnetite  are  present  in  an  oxide 
growing  under  diffusion  control,  not  all  of  the  cations  that  enter 
the  wustite  phase  at  the  Fe/FeO  boundary  contribute  to  the 
growth  of  the  wustite  phase  at  the  FeO/Fe^O^  phase  boundary;  some 
of  the  cations  must  diffuse  further  through  the  magnetite  phase 
to  contribute  to  the  growth  of  the  magnetite  at  the  reaction  front 
(31.  When  this  partitioning  of  cations  between  the  two  phases  is 
included  in  the  analysis,  the  following  relationship  is  obtained 
between  k^tFe^O),  the  parabolic  rate  constant  calculated  with  the 
aid  of  Eq.  (1),  and  kp(FefiO) , the  rate  constant  for  growth  of 
mzgnetite  in  contact  with  a growing  wustite  phase  [3,10]: 


k* (Fe.O) 
P « 


[1  + £V 


Fe^C/ 


'"'F%0)'lVF®60) 


(2) 


where  £ and  n are,  respectively,  the  thicknesses  of  the  growing 
layers  of  Fe^O  and  Fe^O. 

The  equivalent  form  for  Fe^O  is  given  by  [3,10] 


*;<Fevo)  - (l  + «nvHey^v 


(3) 
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and  the  parabolic  rate  constant  for  the  total  oxide  layer 
is  given  by 

k (oxide)  « { [k ' (Fe,  0) ] 1/2  + [k’ (Fex0) ] 1/2}2  (4] 

P p V p 0 

2. 3 Transport  Properties  of  the  Oxides 

The  transport  properties  of  wustite  and  magnetite  that  are 

necessary  to  solve  Eq.  (1)  have  been  redetermined  subsequent  to 

publication  of  ref.  3.  The  tracer  diffusion  coefficients  and  the 

correlation  factors  for  diffusion  of  iron  cations  in  wustite  were 

measured  by  Chen  and  Peterson  [11]  at  approxirately  800,  1000  and 

1200°C.  Data  at  each  temperature  of  interest  in  the  range  of  800 

to  1050°C  were  obtained  [12]  by  interpolating  the  data  of  Chen 

and  Peterson  with  the  aid  of  composition  data  compiled  by  Giddings 

T 

and  Gordon  [5].  The  values  obtained  for  the  quotient  DFe/9Fe» 
which  is  the  self-diffusion  coefficient  for  cations  in  wustite, 
are  listed  in  Table  I. 

The  tracer  diffusion  coefficients  for  iron  cations  in  magnetite 
have  been  measured  as  a function  of  oxygen  partial  pressure  in  the 
temperature  range  of  900  to  1400°C  by  Dieckmann  and  Schmalzried 
[13-15];  their  results  were  extrapolated  to  800  and  850°C.  Magne- 
tite contains  excess  iron  at  low  oxygen  potentials,  but  it  is 
deficient  in  iron  at  high  oxygen  potentials  [13-15].  Diffusion  by 
a vacancy  mechanism  predominates  at  high  oxygen  potentials.  The 
correlation  factor  for  diffusion  by  a vacancy  mechanism  in  magnetite 
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is  gFe  ■ 0.781  (15).  Diffusion  may  take  place  by  an  interstitial 
mechanism  or  by  an  intersticialcy  mechanism  in  magnetite  at  low 
oxygen  potentials  (14).  While  the  exact  mechanism  has  not  yet 
been  determined,  inspection  of  the  spinel  structure  of  magnetite 
suggests  that  in  intersticialcy  mechanism  would  be  energetically 
favorable  relative  to  an  interstitial  mechanism  because  the  latter 
would  involve  very  close  spacings  of  cations. 

If  the  interstitial  cations  in  iron-excess  magnetite  diffuse 
by  an  intersticialcy  mechanism,  then  the  correlation  factor  for 
cation  transport  could  be  less  than  unity.  An  intersticialcy 
irechanisn  could  conceivably  involve  all  combinations  of  normally 
filled  and  empty  tetrahedral  and  octahedral  sites.  One  simple 
possibility,  however,  is  an  intersticialcy  mechanism  that  involves 
only  normally  filled  and  empty  octahedral  sites.  Each  cation  on 
a normal  octahedral  site  is  surrounded  by  six  normally  empty  nearest- 
neighbor  interstitial  sites  that  are  octahedrally  coordinated.  A 
tracer  cation  on  any  one  of  the  six  interstitial  sites  can  displace 
the  cation  on  the  normal  site  to  any  one  of  the  five  empty  inter- 
stitial sites.  Inspection  of  the  structure  suggests,  however, 
that  the  displacement  would  be  collinear.  If  a collinear  inter- 
sticialcy mechanism  is  operative,  the  angle  between  the  jump  vec- 
tors of  the  tracer  cation  and  the  cation  on  the  normal  site  is  zero. 
On  its  next  jump,  then,  the  tracer  has  a greater  than  random  proba- 
bility of  jumping  back  to  its  original  interstitial 
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site,  and  the  motion  of  the  tracer  cation  is  correlated. 


The  corelation  factor  for  a collinear  intersticialcy  mechanism 
for  octahedrally  coordinated  normal  and  interstitial  cations 
in  magnetite  has  not  been  calculated.  Compaan  and  Haven  [16] , 
however,  have  used  an  electrical  analog  technique  to  determine 
correlation  factors  for  a number  of  similar  intersticialcy 
mechanisms.  For  collinear  intersticialcy  diffusion  of  cadmium 
cations  in  CdC^  (rhombohedral  structure) , where  octahedrally 
coordinated  interstitial  cation  sites  have  nearest  neighbor 
cation  sites  in  octahedral  coordination  as  in  magnetite,  both  the 
x-  and  z-compon*;nts  of  the  correlation  factor  are  equal  to  0.4. 

Chen  and  Peterson  [11]  interpret  their  measurements  of  correlation 
factors  in  oxyaen-rich  wustite  in  terms  of  a model  of  Roth  [17] 
and  Koch  and  Cohen  [18]  in  which  clusters  of  vacancies  on  octa- 
hedral cation  sites  group  about  occupied  tetrahedral  sites  in  a 
configuration  resembling  the  structure  of  magnetite.  For  oxygen- 
rich  wustite  at  807  < T < 1200°C  the  measured  g_  AX  5 0.3,  where 
AK  is  the  fraction  of  the  total  translational  kinetic  energy 
at  the  saddle  point  for  a diffusion  jump  [11].  A value  of  0.58  tor 
AK  yields  gfe  ■ 0.78  for  wustite  in  equlibrium  with  iron.  If  AK 
is  assumed  to  be  independent  of  composition,  then  g_  varies  between 
0.5  and  0.59  for  wustite  in  equilibrium  with  magnetite  in  the 
range  of  800  to  1050°C.  Thus,  a low  value  for  the  correlation 
factor  in  iron-rich  magnetite  seems  likely. 

On  the  basis  of  the  above  discussion,  we  assume  that  the 

correlation  factor  for  cation  diffusion  in  magnetite  at  low 

T 

oxygen  pressures  is  0.4.  Thus,  the  expression  for  D_  /g_  in 

r G F6 


-10- 


Eq.  (1)  becomes 


/g 


Fe 


h.  p “2/3  DV  p 2/3 
0.4  1 02  0.781  02 


(5) 


where  the  coefficients  and  Dv  from  Dieckmann  and  Schmalzried 
[15]  are  the  interstitial  and  vacancy  components  of  diffusion  in 
magnetite,  respectively.  Because  gFe  is  less  certain  than  the 
other  parameters  in  the  analysis,  however,  data  are  also  analyzed 
for  comparison  in  terms  of  Eq.  (5)  with  the  factor  0.4  replaced 
by  unity. 

2.4  Molar  Volumes  of  Oxides 

An  average  molar  volume  of  magnetite  was  calculated  using 
a density  of  5.18  x 10  g/m  for  Fe3C>4.  An  average  molar  volume 
of  wustite  at  each  temperature  was  estimated  from  the  data  of 
Hayakawa  et  al.  [19]  for  the  variation  of  the  lattice  parameter 
of  wustite  with  composition  at  950  and  1050°C.  The  average 
molar  volume  of  wustite  was  assumed  to  vary  linearly  with  tempera- 
ture. The  average  values  are  listed  in  Table  I. 


2.5  Transport  Properties  of  Copper 

The  flux  of  oxygen  through  the  product  copper  layer  is  propor- 

Cu  Qu 

tional  to  the  solubility-dif fusivity  product,  NQ  DQ  , where 

Cu 

is  the  atom  fraction  of  oxygen  dissolved  in  copper  and  DQ  is  the 
diffusion  coefficient  for  oxygen  in  copper.  According  to 
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Ramanarayanan  and  Worrell  [20],  the  solubility-dif f usivity  product 
for  the  temperature  range  of  800-1000°C  is  given  by 


„Cu_Cu 
N0  D0 


(8*8-14.,9  x 10”2>VCu  exP(  "-°°-83RTt  8368J)  (6) 


where  VCu,  the  molar  volume  of  copper,  is  equal  to  7.13  x 10_6g/m3. 
The  solubility-dif fusivity  product  determined  by  Pastorek  and  Rapp 
[21]  is  a factor  of  approximately  two  higher  than  that  calculated 

Pn  Cu 

from  Eq.  (6).  Values  of  NQ  DQ  measured  by  Wood,  et  al . [22]  and 
by  Rao  and  Tare  [23]  are  slightly  lower  than  those  calculated  from 
Eq.  (6) . A recent  determination  of  the  maximum  solubility  of 
oxgyen  in  copper  for  temperatures  in  the  range  of  806  to  1050°C  by 
Horrigan  [24]  yielded  values  that  are  70-80%  lower  than  those 
measured  by  Pastorek  and  Rapp  [21].  Thus,  the  solubility-dif fusivity 
product  of  Eq.  (6)  determined  by  Ramanarayanan  and  Worrell  [20] 
seems  to  be  a good  representation  of  the  most  recent  data  available 
in  the  literature.  The  oxygen  partial  pressure  for  equilibrium 
between  copper  and  Cu20  was  obtained  from  the  JANAF  tables  of 
thermodynamic  data  [25] . 

2.6  Parabolic  Rate  Constant  for  Copper 

The  method  employed  to  calculate  the  parabolic  rate  constant 
for  the  growth  of  the  product  copper  layer  was  described  in  ref.  3; 
however,  some  clarification  and  revision  of  the  derivation  is 
required.  The  gross  rate  of  advance  of  the  Cu/CUjO  interface,  x* 
is  given  by  Eq.  (27)  of  ref.  3: 
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2VCuD0UkfU  1/2  F 1/2  kn(Cu) 

C ? S I<2<CU/CU20,  - p"  1/2]  - -V- 


where  kp(Cu)  is  the  parabolic  rate  constant*  for  the  thickening 
of  a layer  of  copper  in  the  hypothetical  case  where  no  copper 
is  left  in  the  aggregate  layer  behind  the  reaction  front  (i.e., 


where  p = 0 in  Fig.  1)  ; 


kg  is  the  Sie verts'  Law  constant. 


The  acutal  rate  of  thickening  of  the  copper  layer  when  iron 
oxide  rods  grow  into  the  copper  layer  is  given  by  Eq.  (28)  cf 
ref.  3: 


at  - x - <3  - « H 


k' (Cu) 


where  kp(Cu)  is  the  parabolic  rate  constant  for  the  thickening 
of  the  copper  layer  when  the  copper  phase  grows  at  the  Cu/Cu20 
interface  but  it  also  exists  in  the  aggregate  region.  The 
actual  rate  constant,  kp(Cu),  depends  on  the  rates  of  growth  of 
all  of  the  product  phases.  The  value  of  k^CCu)  in  Eq.  (7)  depends 
only  on  the  solubility-dif fusivity  product  of  copper  and  on  the 
oxygen  potentials  at  the  boundaries  of  the  copper  layer.  With 
the  insertion  of  x * kp(Cu)/$  and  Eq.  (7),  into  Eq.  (8)  and  use 
of  (de/dt)  * (8/<M  (d«(>/dt)  , Eq.  (27)  in  ref.  3,  one  obtains 


kp' (Cu)  - kp(Cu)/tl  + (1  - f)6/$]  (9) 


*Tms  constant  was  denoted  k'  in  ref.  3,  unlike  the  notation  u»ed 
here  and  in  ref.  10.  " 
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and 


at  u + d - « fi 


k (Cu) 
_E 


♦ 


(10) 


The  value  of  # can  be  calculated  from  Eq.  (33)  of  ref.  3 once  the 

p 

value  of  PI  is  determined. 


3.  RESULTS  AND  DISCUSSIONS 


The  experimental  methods  that  were  employed  to  determine  the 
kinetics  of  the  displacement  reaction  between  iron  and  CU2O  were 
presented  previously  [2,3].  The  thickness  of  the  product  copper 
layer,  the  thickness  of  the  copper  + oxide  product  zone, 

0,  and  the  length  of  magnetite  in  the  oxide  rods,  n»  were  measured 
for  several  reaction  times  at  each  reaction  temperature. 

These  data  were  plotted  as  the  square  of  the  thickness  or 

length  of  a product  phase  vs.  the  reaction  time.  Straight 

lines  were  drawn  through  the  data  using  the  method  of  linear 

2 

regression.  Except  for  the  plots  of  $ vs . reaction  time  at 
800  and  850°C,  the  correlation  coefficients  from  the  linear 
regression  were  generally  greater  than  0.95.  The  straight  lines 
did  not  all  intercept  zero  product  thickness  at  zero  reaction  time, 
which  indicates  that  a single  reaction  mechanism  did  not  prevail 
during  the  entire  course  of  the  reaction.  The  straight  lines 
through  the  data  at  later  reaction  times  demonstrate  that  even- 
tually the  displacement  reactions  were  controlled  by  solid-state 
diffusion  of  the  reactants  through  the  product  phases.  The  slopes 


f 


of  the  straight-line  portions  of  the  plots  are  equal  to  2k' 
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r? 


Values  of  all  the  kinetic  parameters  for  the  Fe/Cu20  displace- 
ment reaction  are  presented  in  Table  II.  The  experimental  values 

of  k' (Fe  0)  were  calculated  from  experimental  values  of  k (oxide) 
p v r p 

and  kp(Fe^O)  with  the  use  of  Eq.  (4).  Experimental  values  of 
(♦/©)  and  (n/C)  were  calculated  using  the  following  equations: 


♦/0  = (kp (Cu)/kp (oxide) ] 


1/2 


(ID 


and 


h/C  * tkp (Fe^O) /kp (Fe^O) ] 


1/2 


(12) 


The  volume  fraction  of  oxide  in  the  three-phase  product  zone  wzs 
calculated  from  Eq.  (24)  of  ref.  3;  i.e.. 


f - (1  + $/e)/(i  + 2vCuVox”1) 


(13) 


:j 


where  VQX  is  a weighted  average  molar  volume  of  the  product  oxide 
phase. 

F 

The  theoretical  values  of  Pn  that  are  presented  in  Table  II 

u2 

are  greater  than  PQ  (Fe0/Fe304)  at  each  temperature  between  800  and 
1050 °C.  Thus,  the  product  oxide  rods  are  expected  to  comprise 
wustite  and  magnetite,  which  was  confirmed  experimentally. 

The  theoretical  and  experimental  values  of  kp^ox*de^  and 
kptFe^O)  agree  well  over  the  entire  range  of  temperature  investigated. 


J 
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i 


i; 

c 
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Note  that  kp(FevO)  is  approximately  equal  to  kp (oxide)  because 
n <<  5 (see  Table  II) . Experimental  and  theoretical  values  of 
kp(FevO)  are  plotted  in  Fig.  2.  The  activation  energies  obtained 
from  the  slopes  of  the  plots  of  log  kp(FevO)  vs.  T are  157,400 
and  166,500  J/mole  for  the  experimental  and  theoretical  rate  con- 
stants, respectively. 

Experimental  and  theoretical  rate  constants  for  the  growth 

of  the  product  copper  layer  agree  well  within  the  temperature 

range  of  approximately  950-1050°C.  Below  ~950°C  the  experimental 

and  theoretical  values  of  kp(Cu)  diverge , the  divergence  increasing 

as  the  temperature  decreases  as  illustrated  in  Fig.  3.  The  straight 

line  through  the  data  plotted  in  Fig.  3 was  determined  by  linear 

regression;  the  activitation  energy  determined  from  the  slope  of 

the  line  is  315,300  J/mole.  The  theoretical  activation  energy 

is  240,500  J/mole.  The  value  of  k' (Cu)  depends  both  on  the  rate 

P 

at  which  oxygen  diffuses  through  copper  and  on  the  rate  at  which 
the  product  oxide  rods  grow  into  the  product  copper  layer  [see  Eq. 
(8)].  The  activation  energy  from  the  slope  of  a plot  of  log  kp(Cu) 
vs.  T * is  an  activation  energy  for  the  net  process  of  growth  of  the 
copper  layer,  and  it  should  not  be  compared  solely  to  the  activation 
energy  for  the  transport  of  oxygen  in  copper  as  obtained  from  Eq. 

(5). 

The  temperature  dependencies  of  the  experimental  and  theoretical 
rates  of  growth  of  the  product  oxide  rods  agree  fairly  well.  Thus, 
the  difference  in  the  temperature  dependencies  of  the  experimental 
and  theoretical  values  of  k^(Cu)  appears  to  be  related  mainly  to  the 
dependence  on  temperature  of  the  flux  of  oxygen  through  copper. 


J 
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Better  agreement  between  experiment  and  theory  is  obtained  if  the 

absolute  magnitudes  of  the  exponential  and  preexponential  terms  in 

the  expression  for  the  solubility-dif fusivity  product  for  oxygen 

in  copper,  Eq.  (6)  , are  increased.  Using  the  appropriate  limits 

of  precision  in  Eq.  (5)  one  obtains  an  . activation  energy  of 

255,000  J/mole  for  the  theoretical  value  of  k'(Cu).  The  rate 

P 

constant  is  increased  slightly  at  1050°C  and  decreased  slightly 

at  800°C.  The  small  change  in  the  solubility-dif fusivity  product 

has  no  effect  on  P_  or  on  the  other  kinetic  parameters.  In  view 

°2 

of  the  high  sensitivity  of  kp(Cu)  to  the  transport  properties 

of  the  product  phases,  the  agreement  between  experimental  and 

theoretical  values  of  k' (Cu)  is  considered  to  be  good. 

P 

The  experimental  and  theoretical  values  of  $/8  are  generally 
in  good  agreement.  The  variations  in  the  experimental  and  theoreti- 
cal values  of  $/0  with  temperature  follow  approximately  the  same 
trends  as  k^(Cu)  and  kp (oxide)  because  of  the  dependence  of  <j>/0 
on  these  rate  constants  [Eq.  (11)]. 

The  experimental  and  theoretical  values  of  f,  which  were 
calculated  from  Eq.  (13) , are  in  good  agreement  at  all  temperatures 
in  the  range  of  800  to  1050°C.  Note,  however,  that  f is  approximately 
equal  to  (1  ♦ 2VCuvFe  0)-1  when  n <<  5 and  4>  <<  0,  which  is  found 
experimentally.  Therefore,  the  good  agreement  between  experimental 
and  theoretical  values  of  f simply  means  that  the  theory  is  able  to 
predict  n <<  £ and  $ <<  0. 


-17- 


Por  temperatures  between  950  and  1050°C,  the  theoretical  values 
of  kp(FegO)  are  in  excellent  agreement  with  the  experimental 
values  of  kptFe^O).  The  experimental  and  theoretical  values  of 
the  rate  constant  differ  sharply  for  temperatures  below  950°C, 
as  shown  in  Fig.  2.  The  theoretical  values  of  n/S  agree  well  with 
the  experimental  values  of  n/5  for  temperatures  between  950  and 
1050 #C,  but  the  two  values  differ  markedly  for  T ? 950°C,  as 
shown  in  Fig.  4.  The  variation  of  n/t  with  T is  simply  a 
reflection  of  the  temperature  dependencies  of  k^CFe^O)  and 
k^(Fev0)  because  of  the  dependence  of  n/C  on  these  rate  constants 
according  to  Eq.  (12) . 

p 

The  theoretical  values  of  Pft  for  temperatures  above  950°C  are 

°2 

less  than  the  oxygen  partial  pressures  at  which  the  mechanism 

of  diffusion  in  magnetite  changes  from  one  involving  primarily 

\ 

interstitial  cations  to  one  involving  primarily  vacant  cation  sites. 
Thus,  the  growth  of  magnetite  during  a displacement  reaction 
between  iron  and  Cu20  above  950°C  -is  supposed  to  occur  by  dif- 
fusion of  iron  cations  through  magnetite  by  an  intersticialcy 
mechanism. 

The  use  of  a correlation  factor  of  unity  for  an  interstitial 

or  intersticialcy  mechanism  yields  theoretical  values  of  kp(Fe^O) 

that  are  seven  to  eight  times  lower  than  the  experimental  rate 

constants  at  temperatures  above  950°C  (see  Fig.  2) . The  values 
F 

of  P_  in  the  latter  case  are  only  a factor  or  two  lower  than 
°2 

in  the  case  in  which  the  correlation  factor  is  taken  as  0.4.  Thus, 
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the  generally  good  agreement  bewteen  the  other  experimental  and 
theoretical  kinetic  parameters,  except  n/C  » is  maintained  when 
the  correlation  factor  is  taken  as  unity. 

For  temperatures  below  950*C,  the  general  mechanism  of  the 
displacement  reaction  between  iron  and  Cu^O  is  qualitatively 
the  same  as  at  higher  temperatures;  however,  the  disagreement 
between  theory  and  experiment  for  magnetite  growth  at  these  lower 
temperatures  indicates  that  the  mechanism  for  growth  must  differ. 
Several  reasons  could  account  for  the  change  in  the  mechanism  of 
growth  of  magnetite  at  lower  temperatures.  For  example,  a change 
from  predominant  lattice  diffusion  of  cations  through  magnetite 
at  high  temperatures  to  predominant  interfacial  diffusion  of  cations 
at  lower  temperatures  could  yield  the  observed  dependence  of 
kp(Fe^O)  on  temperature.  Also,  the  lower  activation  energy  obser- 
ved at  temperatures  below  950°C  is  characteristic  of  the  activation 
energy  for  the  growth  of  magnetite  when  migration  of  cations 

through  magnetite  by  a vacancy  mechanism  makes  a major  contribution 

F 

to  the  growth  process.  Values  of  P that  are  greater  than  the 

°2 

oxygen  partial  pressure  for  the  change  from  an  intersticialcy  to 
a vacancy  mechanism  of  diffusion  in  magnetite  could  produce  the 
desired  effect.  For  example,  values  of  PQ  necessary  to  yield 
agreement  between  experimental  and  theoretical  values  of  kp(Fe^O) 
at  temperatures  below  950°  are  very  close  to  the  partial  pressures 
of  oxygen  for  equilibrium  between  copper  and  CUjO.  The  required 

p 

values  of  P_  are,  however,  greater  than  the  oxygen  potentials 
u2 
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necessary  to  form  CuFe02  and  they  would  yield  very  poor  agreement 
between  the  experimental  and  theoretical  values  of  kp(Cu).  It 
is  noted  that  the  morphology  of  the  product  zone  of  the  Fe/Cu20 
displacement  reaction  at  low  temperatures  does  not  conform 
exactly  to  the  model,  i.e.,  at  low  temperatures  the  thickness  of 
the  product  copper  layer  is  not  much  greater  than  the  spacing 
between  the  product  oxide  rods,  as  is  illustrated  in  Fig.  5.  At 
low  temperatures,  the  assumption  of  unidirectional  diffusion  of 
oxygen  atoms  through  the  product  copper  layer  is  apparently  not 
valid.  Thus,  the  kinetic  model  cannot  be  applied  unambiguously  to 
the  Fe/CUjO  displacement  reaction  at  low  temperatures.  Which,  if 
any,  of  these  reasons  account  for  the  change  in  the  mechanism  of 
the  growth  of  the  magentite  phase  cannot  be  decided  at  present. 


4 , CONCLUSIONS 

The  products  of  the  displacement  reaction  between  iron  and  Cu20 
for  reaction  temperatures  between  800  and  1050°C  are  copper, 
wustite  and  magnetite.  Excellent  agreement  is  obtained  between 
experimental  and  theoretical  kinetic  parameters  for  reaction  tem- 
peratures between  950  and  1050°C  if  the  mechanism  of  diffusion 
of  iron  cations  in  magnetite  for  displacement  reactions  above  950°C 
is  an  intersticialcy  mechanism  with  a correlation  factor  of  0.4. 

The  mechanism  of  the  growth  of  magnetite  changes  at  approximately 
950°C.  The  mechanism  of  the  displacement  reaction  between  Fe 
and  CUjO  at  lower  temperatures  has  not  been  determined  although 
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interfacial  diffusion  of  cations  and  enhanced  lattice  diffusion 
of  cations  by  a vacancy  mechanism  in  magnetite  at  lower  tempera- 
tures may  be  important  factors . 
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QUANTITIES  NECESSARY  TO  CALCULATE  THE 


RATE  OF  GROWTH  OF  WUSTITE 


C 

V 

K 

(m3/mole)xl05 

(DFe/9Fe) -m2/s 

800 

0.9227 

1.230 

7.31  x 10“15Pn  ~0*118 
u2 

850 

0.9194 

1.233 

1.28  x lO-13?^  ~0-070 
U2 

900 

0.9165 

1.236 

1.34  x 10"12Po  “°*027 

950 

0.9137 

1.239 

9.30  x 10"12Po  °*011 
°2 

1000 

0.9111 

1.242 

4.55  x 10"11?-.  °*044 
°2 

1050 

0.9088 

1.245 

1.80  X 10"10P.  °-076 
°2 

TAliLi;  cAi  tiii.-a.uJAL  AND  fHEURETlCAI.  KINETIC  PARAMETERS  FOR  THE  Fe/Cu,0  DISPLACEMENT  REACTION 
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FIGURE  1.  A model  for  the  displacement  reaction  between  iron 
and  Cu20.  The  relative  dimensions  of  the  product 
phases  are  greatly  exaggerated  for  the  purpose  of 
illustration.  It  is  assumed  that  <p  >>  (t  + o) 
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2.  The  dependence  on  reaction  temperature  of  the  experi- 
mental and  theoretical  parabolic  rate  constants  for  the 
growth  of  wustite  and  magnetite. 
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FIGURE  3.  The  dependence  on  reaction  temperature  of  the 
experimental  and  theoretical  parabolic  rate 
constants  for  the  growth  of  the  product  copper 
layer . 
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FIGURE  4.  The  dependence  on  temperature  of  the 
experimental  and  theoretical  values 
of  (n/e). 
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